CMB Constraint on Radion Evolution in the Brane World Scenario 



K. C. ChaiE 

Department of Physics and Institute of Theoretical Physics, 
the Chinese University of Hong Kong, Hong Kong SAR, People's Republic of China and 
Department of Physics, New York University, 4 Washington Place, New York 10003, USA 

M.-C. Chi]!! 

Department of Physics and Institute of Theoretical Physics, 
the Chinese University of Hong Kong, Hong Kong SAR, People's Republic of China 

(Dated: February 5, 2008) 

In many versions of brane model, the modulus field of extra dimensions, the radion, could have cos- 
mological evolution, which induces variation of the Higgs vacuum expectation value, {H), resulting 
in cosmological variation of the electron mass me. The formation of Cosmic Microwave Background 
(CMB) anisotropies is thus affected, causing changes both in the peaks positions and amplitudes 
in the CMB power spectra. Using the three-year Wilkinson Microwave Anisotropy Probe (WMAP) 
CMB data, with the Hubble parameter Hq fixed to be the Hubble Space Telescope (HST) result, 
72 km s"^ Mpc"\ we obtain a constraint on p, the ratio of the value of {H) at CMB recombination 
to its present value, to be [0.97, 1.02]. 

PACS numbers: 06.20.Jr, 98.80.Cq, 98.70.Vc 



I. INTRODUCTION 

Investigations in string/M theory suggest the existence 
of extra dimensions. In particular, we may live in a brane 
world in a bulk of dimensions greater than 4. In this 
brane world picture, the matter fields are confined on the 
brane which is a (1-f 3)-dimensional hypersurface embed- 
ded in the bulk, while gravity can propagate freely in 
the whole extended universe. For a review on extra di- 
mensions and brane worlds see e.g. [l|, By restricting 
the standard model particles on the brane, brane world 
models have evaded the strong constraints on the size of 
extra dimensions. Up till now Newtonion gravity is only 
tested to submillimeter scale [1], and existence of extra 
dimensions below this scale remains possible. From the 
phenomenological point of view, extra dimensions open 
a new door to address the hierarchy problem, the large 
scale discrepaiicy between the Planck scale and the elec- 
troweak scale [J]. For example, in the Randall-Sundrum 
two-brane model Q , one can obtain a large scale discrep- 
ancy by fine-tuning the interbrane distance, the radion; 
thus some mechanisms have to be invoked to stabilize the 
radion [1, 0|- 

There are many moduli fields in string theory and the 
radion is an example. The moduli fields have to be sta- 
bilized so as not to upset the limits from long range force 
experiments and tests of general relativity. Other than 
invoking stabilization potentials, there could also be some 
attractor mechanism that drives the theory towards gen- 
eral relativity ^ . The radion may also act as a chameleon 
field and acquire mass by self-interactions to avoid the 



'Electronic address: kcc274@n yu.edii] 

t Electronic address: ,mcchu@phy.cuhk.edu.hk] 



most restrictive of current bounds [1, [T3| • If these sce- 
narios hold, it is possible that the radion still evolves 
after the epoch of Big Bang Nucleosynthesis (BBN) . 

Studies of brane world models can be divided into two 
main streams. On one hand, more complex models are 
studied in order to strengthen the connection between 
string theories and brane models; on the other hand, 
efforts are made to confront brane models with obser- 
vational data pd| . In the brane world scenario, evolv- 
ing radion induces variation of Higgs VEV, and therefore 
constraining the evolution of the radion implies limiting 
the evolution of Higgs VEV. In Ref. [H, [ll, the con- 
straints on the evolution of the radion are obtained from 
BBN and the best fit to the observational abundances of 
D, "'He and ''Li suggests a small variation in the radion 
during BBN compared to its present value. Molecular 
spectral lines can be used to determine fj, = mp/nie p^ . 
Recently, Reinhold et al. used quasar H2 spectral lines to 
obtain A/x//x = (2.4 ± 0.6) x 10"'^ for a weighted fit 
The constraints on /i can be interpreted as an indication 
that the Higgs VEV (and hence the radion) at redshift 
2-3, was different from its value today. The Cosmic Mi- 
crowave Background (CMB) anisotropies are measured 
with high precision by Wilkinson Microwave Anisotropy 
Probe (WMAP) [l^ and other observations. Thus it is 
timely to study the effects of the existence of extra di- 
mensions on CMB and fill the "redshift gap" between 
BBN and quasar constraints using data from CMB. In 
this article we shall use the three-year WMAP data to 
constrain the evolution of the radion. 

Within the brane world scenario, an extra term, called 
the dark radiation term, appears in the modified Fried- 
mann equations [it'I , and it causes great difficulty in cal- 
culating the Sachs- Wolfe Effect j^,i9„M, 21]. Nonethe- 
less, in this work we will constrain the radion evolution 
through the variation of constants caused by the vari- 
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ation of Higgs VEV ignoring the dark radiation term. 
CMB has been used to constrain the Higgs VEV using the 
pre-WMAP data in [23|- However it is desirable to get 
an improved bound using the high precision three-year 
WMAP data, including the polarization power spectrum. 
We also demonstrate that the constraints can be tight- 
ened substantially if we make use of HST measurement 
of Ho. 

This rest of the paper is organized as follows. In Sec- 
tion [Hi we carry out the conventional dimensional reduc- 
tion to get the radion dependence of various constants. 
We shall see that only the Higgs VEV is expected to vary 
in the context of brane world models. We then describe 
the effect of variation of Higgs VEV on the CMB power 
spectra in Section IHI Al The numerical constraints on 
the radion (or Higgs VEV) are presented in Section llllBI 
We summarize in Section HVl 



II. DIMENSIONAL REDUCTION AND LOW 
ENERGY EFFECTIVE ACTIONS 

In this section we will start from the higher dimen- 
sional action and carry out the dimensional reduction to 
derive the dependence of the constants on the radion. 
There are many works on effective actions in extra di- 
mensions, e.g. [1, [l^, [23 | . a nd the formalism here fol- 
lows closely that in Ref. [l^, [2^. We shall dimension- 
ally reduce the higher dimensional gravitational action 
to obtain its low energy effective action first. The higher 
dimensional action reads 
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Ricci scalar ^^^R, S can be written as 



S = 
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where we have defined 
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with Vq being the volume of the extra dimensions today. 
We shall work in the Einstein frame with pure Ricci scalar 
in the gravitational action. To do so we will apply the 
conformal transformation 
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with 
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For concreteness, we further assume that the extra di- 
mensional manifold is homogeneous and isotropic. Hence 
the extra dimension metric takes the simple form: 



diag(fe2,62,...,52). 



(7) 



After the transformation Eq. [5] and using the ansatz 
Eq. [7l the action S reduces to 



S = J d^x^—g 



1 n(n + 2) 1 



(8) 



To make the scalar field canonical, we define a new scalar 
field, the radion a, as 



where n denotes the number of compact extra dimen- 
sions, and represents the bulk spacetime coordinate, 
A = 0,l,...,3-|-n. The higher dimensional Ricci scalar 
is denoted by i?, and G is the determinant of the 
full spacetime metric Gab- C,m is the matter field La- 
grangian density, which may include scalar fields, vector 
boson fields and Dirac fermion fields. 

To proceed we shall take the metric ansatz 

ds^ = GABdX^dX^ = g^^{x)dx^'dx'' + h,j{x)dy'dy^ . 

(2) 

In this metric the Greek indices run over 0, 1, 2 and 3, 
while the Latin indices i and j run over the extra di- 
mensions from 4 to 3 -|- n. We are only interested in 
the zero-mode of the Kaluza Klein expansion, and so the 
extra dimensional metric hij{x) does not depend on the 
extra dimension coordinate. Furthermore, we have as- 
sumed that the metric in Eq. [2] is block-diagonal because 
the vector-like conection GJ^ vanishes for zero-mode [2^. 
The extra dimensions are compactified on an orbifold and 
the dimensionless coordinate assumes values in the in- 
terval [0,1]. Expressing (4-i-n)/j [j^ ^]-^g terms of the 4D 
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Finally the effective action reads 
S = I d^xj- 



(10) 

We now study the matter sector in more details. Al- 
though the matter fields only live on the brane, their 
actions are still affected by the existence of extra dimen- 
sions because of the conformal transformation Eq. [5] Let 
us begin with the minimally coupled scalar field 
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After the transformation Eq. [51 we have 

1 
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where we have defined 

$ — exp —^4(7 



2(n + 1) ' 

to make the kinetic term canonical. 

If 1^(0) is taken to be a simple renormalizable potential 
of the form 



(14) 



we then have 



V{^) = \^^^^n^as^^^^^^'' ^\v'^\ (15) 

We see that the mass of the scalar field is dependent 
on the radion in the Einstein frame for both signs of 
\j? . One of the most important scalar fields in particle 
physics is the Higgs field. Masses of the fermions and 
quarks can be generated by their Yukawa couplings to 
the Higgs field. In this mechanism, the fermion masses 
are proportional to the Higgs VEV, (if). Thus we get 
the radion dependence of the fermion mass 
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For the gauge field, the action is given by 
1 
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where F^^, is the gauge- invariant field strength tensor. 
However the action is invariant under the conformal 
transformation Eq. [5l and so the effective 4D coupling 
constant is the same as the higher dimensional g^. 

Similar techniques can be applied to the Dirac field tp 
with mass m 1261: 



5'Dirac = J d^xy/^ (iVe'^7»-D^V' " wV'V') , (18) 

where e*'' is the vierbein and is the covariant deriva- 
tive. Applying Eq. [5] and redefining the field tp as 



= exp 

we get the canonical action 
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We can read out the radion dependence of the fermion 
mass, which is the same as Eq. 1161 

Thus, for the brane world models in the Einstein frame, 
only the fermion masses, among other fundamental con- 
stants, acquire radion dependence and are expected to 
vary. Or, in the framework of Standard Model, we may 
say that only (H) is radion-dependent while the Yukawa 
coupling is constant. 



III. NUMERICAL CONSTRAINTS ON THE 
EVOLUTION OF THE RADION 

In the last section, we see that only (H) is expected to 
vary in the brane world scenario. Here we first discuss 
the effects of variation of (H) on the CMB power spectra, 
and then we present the numerical constraints on (H) 
and the radion using the the three-year WMAP data. 



A. EfTects of variation of Higgs VEV on CMB 
power spectra 

The variation of (H) induces changes in the Fermi con- 
stant Gf, the quark masses, the nucleon binding energy 
and the electron mass nie [23] ■ The variation oi Gp is 
not relevant in the epoch of CMB recombination. We 
can also ignore the changes in strong nuclear force and 
nuclear binding energy caused by the variation of quark 
mass. Since the baryon mass is dominated by the QCD 
scale parameter Aqcd, and the quark masses contribute 
at only a few percents level, we shall neglect the inffuence 
of (H) on the baryon mass. Thus the variation in (H) 
boils down to variation of me [28j . which modifies the 
recombination history. 

The variation of rUe enters CMB through the bind- 
ing energies of hydrogen. Hi, and helium. Hen, Thom- 
son cross-section ctt, and the recombination coefficients 
a and the two-photon decay rates of hydrogen and he- 
lium. In the Appendix, we list the evolution equations 
and modifications to take care of variations of mg. Here 
we only sketch the modifications to be made. The bind- 
ing energies scale as mg, ctt is proportional to m~^, and 
the two-photon decay rates vary as nie [2§|. For a, one 
can derive a differential equation relating the matter tem- 
perature Tm and nie [2^ , through which the dependence 
of a on TTie can be deduced using the empirical fitting of 
a as a function of Tm in the literature. 

Among these effects caused by the variation of nie, the 
change in the binding energy of hydrogen is most signifi- 
cant on the power spectra; less important is the effect of 
Thomson cross-section, while the effects of the recombi- 
nation coefficients and two-photon rates are small. 

For convenience, we define 



(H) 
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where (-ff)cMB and {H)q denote the Higgs VEV in the 
era of CMB recombination and today respectively. 

Shown in Fig. [T] is the CMB temperature power spec- 
trum with p being 1, 1.05 and 0.95 respectively with the 
other parameters being the standard ones. We see that 
both the positions and amplitudes of the peaks change 
for p = 1.05 and 0.95 in comparison with p = 1. For 
p — 1.05 (0.95), because the binding energy cx nie and it 
dictates the epoch of CMB recombination, CMB recom- 
bination takes place at an earlier (later) time. Therefore 
the sound horizon is smaller (larger), and the distance to 
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FIG. 1: The CMB temperature power spectrum with p—1, 
1.05 and 0.95 respectively. Other cosmological parameters 
assume the standard values. 



FIG. 2: The CMB _E-polarization power spectrum with p= 
1.05 and 0.95 respectively. 



the last scattering surface is larger (smaller); the peaks 
shift to the larger (smaller) I scales. The change in ampli- 
tudes is due to early Integrated Sachs- Wolfe (ISW) effect 
and damping effect. For p = 1.05 (0.95), the residual ra- 
diation is higher (lower) right after recombination, which 
gives an enhanced (reduced) early ISW effect and hence 
the first peak is boosted (diminished). Early recombina- 
tion also means greater Hubble rate and narrower visi- 
bility function at decoupling time, and thus the damping 
is reduced and the power in high I scales is enhanced. 

Shown in Fig. [2] is the i?-polarization power spectrum 
with p—\, 1.05 and 0.95 respectively. The main fea- 
tures of this power spectrum can be understood with the 
tight coupling semi-analytic formula for the polarization 
strength Op; of a particular k mode (soj : 



Qpi{k) 



5/c6i(fc, 77*) 



6f(r/,) [A:(?7o - ??*)] 



,]i{k{r]^-V*)), (22) 



where 770 and 77* are the conformal times of today and 
the last-scattering surface respectively. &i{k,'ri^) is the 
dipole of the photon distribution. Again, for larger 
r/igj recombination occurs at an earlier conformal time, 
and because the spherical Bessel function peaks at I ~ 
k{rio — rj^,), the peaks shift to the larger Is. The power of 
the peaks is higher because the dipole, which is like the 
velocity of a fluid in the tight coupling limit, is larger at 
the earlier time. 

In the above discussions, we have assumed that rue 
is constant in the whole calculation. Because the effect 
of TTie on CMB is only important near the decoupling 
time, even if we assume that rrie evolves, e.g. as a linear 
function of the scale factor a, only the value of tti^ at the 
decoupling time matters. Thus CMB indeed gives the 
constraint on to. at 2 ~ 1000. 



B. Numerical constraints on the Higgs VEV and 
radion by CMB 

In this section, we constrain the range of (H) (or ra- 
dion) using the three-year WMAP data To do 




FIG. 3: The marginal distributions of the free parameters, 
constrained by the three-year WMAP data. In addition to 



the free parameters Wb = 0,bh , ujc ~ i^ch , Hq, Zre 



A, 



and p, shown also are the derived distributions of the density 
parameter of matter {Q,m) and cosmological constant {Q,a). 
The Hubble parameter is allowed to vary in the range from 
40 to 100 km Mpc~^ (we will suppress this unit after- 
wards). Here and thereafter, the maxima of the distributions 
are arbitrarily normahzed to 1. 



SO, we make use of the Markov Chain Monte Carlo 
(MCMC) method implemented by the engine CosmoMC 
[3l| . which searches for the maximum of the likelihood 
function. The theoretical CMB spectra are calculated by 
the Boltzmann code CMBFAST 32]. We vary the fol- 
lowing set of parameters: the Hubble parameter, Hq, the 
baryon density, ujb — Vt^h? {h — Hq/XQQ km s^^ Mpc~^), 
the cold dark matter density, ujc — ^ch^, the reioniza- 
tion redshift, Zio, the primordial fluctuation amplitude. 
As, the spectral index, n^, and p. Flatness of the universe 
is assumed in all the calculations. 

The marginal distributions of the parameters are 
shown in Fig. [3] In particular, the 95% confidence inter- 
val (C. I.) for p is [0.85, 1.11]. The constraint is relatively 
weak. Most of the other standard cosmological parame- 
ters are well constrained in the usual ranges, except the 
Hubble parameter. Although the mean of Hq being 67.6 
is still close to the "canonical" value 72, its distribution 
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FIG. 4: The contour marginal distributions of cosmological 
parameters plotted against p. Among the free parameters, 
Ho is the most strongly degenerate with p. 



FIG. 6: The derived marginal distributions of Inp = 
-^■\/^(o"CMB - o-o) for 40 < Ho < 100 (solid line) and 
Ho = 72 (dashed line) respectively. Note that the latter has 
much smaller spread. 



0.02 0,021 0.022 0.023 0.024 O.OB 0.1 0.12 




FIG. 5: Same as Fig. |3l but with Ho fixed to be 72. 
spread of p is much smaller than that in Fig. |31 



The 



spreads wide. This suggests that Hq may be degener- 
ate with p. To see the degeneracy between p and other 
parameters explicitly, we display the contour marginal 
distributions of various cosmological parameters plotted 
versus p in Fig. [H Among the parameters, ujb and ujc 
show slight degeneracy with p, but p is strongly degener- 
ate with Hq. This suggests that we may obtain a tighter 
bound if we use the measurement of by the Hubble 
Space Telescope (HST) Key Project to break the degen- 
eracy between and p. We carry out the MCMC run 
once again but with = 72 [33!|. Indeed, the bounds 
on p are tightened substantially, as can be seen in Fig. [S] 
The 95% C. I. now becomes [0.97, 1.02]. 

Using Eq. [111 we get the change of the radion 



fCMB — ctq 



2n 
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Inp, 
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which is proportional to In p. In Fig. [SI we show the de- 
rived distribution of In p. One sees again that the spread 
of the distribution of (ctcmb — fo) for Hq being allowed 
to vary from 40 to 100 is much larger than that for fixed 
Hq = 72, and their 95% C. I. are [-0.094, 0.17] and [- 
0.019, 0.034] respectively. 



The constraint on Higgs VEV by CMB is first deduced 
in [l^] using pre-WMAP data sets. There is no marked 
difference between their constraint 0.92 < p < 1.13 
and ours with Hq being allowed to vary in the interval 
40 < i/o < 100. The small difference may be due to the 
use of different data sets and/or statistical methods. The 
data set we used is the recent three-year WMAP data 
with the polarization spectrum [l^ while the pre-WMAP 
temperature power spectrum is used in [12] ■ They re- 
stricted -ffo in the interval [50, 80] and more importantly 
VLm in the interval [0.3, 0.4], which is a small portion of 
the distribution in Fig. [SI We also allow for two more 
free parameters, Zre and As in our MCMC runs. That's 
why we get bounds comparable to theirs although we 
use more accurate data set. Moreover, we show that the 
bounds on p can be tightened substantially if the HST 
measurement of is used. 

Since the upcoming Planck satellite mission is going to 
measure the temperature power spectrum to as high as 
I ^ 2500 and the £'-polarization spectrum to / ^ 1500, we 
expect that there will be tremendous improvement in the 
constraint on p. We can forecast the improvement that 
Planck will bring quantitatively using the Fisher matrix, 
which has been widely used to predict the expected un- 
certainties in future experiments (see e.g. |34|). Under 
the assumption of Gaussian perturbations and Gaussian 
noise, the Fisher matrix takes the form 



; x,Y 



dpi 



dpj 



(24) 



where pi is the ith free parameter and Cxi is the Ith 
multipole of the observed spectrum of type X, which 
can be the temperature, temperature-polarization and 
E-polarization spectra. The experimental precision is 
encoded in the covariant matrix Covixy- With the ex- 
pected Planck power spectra, the constraint on p is tight- 
ened by a factor of 7 when Hq is allowed to vary; the 
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current constraint with Hq fixed will be tightened by a 
factor of 5. 

We compare the present CMB constraints with other 
constraints in Table |T1 We note that the BBN constraint 
and the quasar absorption line constraint are of opposite 
signs. Even if both constraints are found to be correct, 
that does not immediately rule out the evolution of the 
radion since the radion may stabilize to its present value 
in an oscillatory manner. If this is the case, we may 
find that the radion takes a value very close to 1 near 
the epoch of CMB recombination. Unfortunately, the 
current CMB data is not discriminating enough. With 
the forthcoming Planck data, the CMB constraint will be 
tightened substantially and may tell us if this interesting 
scenario holds or not. 



TABLE I: The present constraints on the variation of {H) at 
various redshifts. The CMB constraint fills the gap in the 
"redshift ladder" in between BBN and quasar. The future 
Planck data will tighten the CMB constraint ( so that it is 
compatible with other constraints.) 



Redshift z 


Observations 


95 % C. I. for {{H), - 
{H)o)/{H)o, where {H), de- 
notes the value of Higgs VEV 
at redshift z 


10"' 


BBN [Hi 


[0.00,0.04] 


1000 


CMB 


Ho free: [-0.15,0.11]; Ho=72: 
[-0.03, 0.02] 


2-3 


Quasar absorp- 
tion lines 15] 


Weighted fit: -0.9 x [1.8 x 
10"^3.0 X 10"^] " 



"The prcfactor 0.9 is due to our assumption that 10% of the proton 
mass is contributed by the rest mass of the quarks. 



IV. CONCLUSION 

In the brane world scenario, the evolution of radion 
induces variation of (H). Among its consequences, we 
expect only the variation of rrie to be relevant in the 
epoch of CMB recombination. Variation of me changes 
the time of CMB recombination, causing changes in both 
the peak positions and amplitudes of the CMB power 
spectra. Thus we can constrain the radion evolution via 
limiting the variation of (H) using the CMB data. With 
the three-year WMAP data and Hq as a free parameter, 
we obtain the 95% C. I. for p to be [0.85, 1.11]; when 
we fix Ho to be the HST result 72 km s"^ Mpc'^ the 
constraint for p is tightened to be [0.97, 1.02]. In terms 
of In p, which is proportional to the change in radion a by 
Eq. El the corresponding 95% C. I.'s are [-0.094, 0.17] 
and [-0.019, 0.034] respectively. Although the current 
CMB data is not discriminating enough, the upcoming 
Planck data should tighten the bounds by a factor of 5 
or so. 



Appendix 

In this appendix we list the evolution equations for cal- 
culating the ionization history with the effect of modified 
nie taken into account. The detailed evolution of the ion- 
ization fraction Xp can be modeled by the simple Peebles 
recombination [3 51 ] , which is basically a two-level approx- 
imation, or using RECFAST [36|] , the evolution equations 
of which have been checked against more detailed multi- 
level calculations. We use a modified RECFAST which 
allows for variation of me to evaluate the evolution of Xe, 
as in Ref. [13] • The ionization history can be modeled by 
evolving the ionization fraction of H, Xp, ionization frac- 
tion of Hci, xhcii 
[3^ for details): 



and the matter temperature Tm (see 



dz 



H{z)a + z) 

/3h(1 - a:p)exp ( - 



[xeXpHaail 



kT, 



M 



with 



Xp j 



(25) 



(26) 



dx 



Hot 



Cneii 



dz 



a;Hciia;e'^HaHci 



H{z){l + z) 
/3hci(/hc - a^Hcii) exp 



.(27) 



with 

Chcii 

and 

dTM 
dz 



1 -I- ii'HeiAHe"H(/Ho " a^Hon) exp ^- 



1 + K]ici{Anc + /3hoi)"-h(/ho - shch) cxp ( - '"'"fcr'^'^ 



(28) 



3'mecH{z){l + z) I + /hc + X, 



■{Tm-Tu 



2T, 



M 



l + Z 

(29) 

The meanings of the symbols are in order. In these 
equations, z is the redshift and H{z) is the Hubble ex- 
pansion rate at z, and hh is the total hydrogen number 
density. Tr denotes the radiation temperature To{l + z). 

i^ii2s — c/121.5682 nm is the Lyman a fre- 
quency. For helium, j^hci2s = c/60.1404 nm, and 
^'Hci2p-2s denotes c/58.4334 nm - c/60.1404 
the frequency difference between Hei 2^p and 



nm, 
Hei 



(121.5682 nm)V[87riJ(z)] and Kho 



(58.4334 nm}^ / [8tt H (z)] give the amount of cosmologi- 
cal redshiftings of 2^p to l^s photons. These frequency 
V dependent quantities need to be modified as v ex. rUe- 

Ah and Ahc denote the two-photon decay rates from 
2^s to l^s for H and He respectively, and they scale as 
me [H. 
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The Thomson cross-section ax scales as ^ as ctt = 

aH (ofHci) is the recombination coefficient of Hi (Hen), 
and the photo-ionization /3h (/^Hou) is related to it by 



= a 



(3 = a 



3/2 



exp 



The recombination coefhcient a can be expressed as 



(30) 



In the literature a is usually parametrized as a function 
of T]\j in the fitting formula. The me dependence of a 
can be extracted from the fitting formula using Eq. [32] 
and hence also the rrie dependence of p. 



a = ^87r(2?+l) 



n.l 



kT, 



M 



2nmf 



exp 



/ B„ 



V ^Tm J J ey — I 



(31) 

where Bn is the binding energy of the nth state and ani is 
the ionization cross-section. The asterisk indicates that 
the sum should be regulated. Using the fact that _B„ 
scales with and the ionization cross-section scales as 
n^r^ [III, we can derive [111 



(32) 
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